Macroscopic quantum tunneling and phase diffusion in a La 2 _ x Sr :r Cu04 intrinsic 

Josephson junction stack 
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We performed measurements of switching current distribution in a submicron I^-rzSrzCuCU 
(LSCO) intrinsic Josephson junction (IJJ) stack in a wide temperature range. The escape rate 
saturates below approximately 2 K, indicating that the escape event is dominated by a macroscopic 
quantum tunneling (MQT) process with a crossover temperature T* w 2 K. We applied the theory 
of MQT for IJJ stacks, taking into account dissipation and the phase re-trapping effect in the 
LSCO IJJ stack. The theory is in good agreement with the experiment both in the MQT and in 
the thermal activation regimes; thus it is confirmed that the enhancement of MQT rate due to 
interaction between junctions in the stack is a universal phenomenon in IJJs. 

PACS numbers: 74.50. +r, 85.25.Cp, 74.81. Fa, 74.72.-h 



I. INTRODUCTION 

Stacks of intrinsic Josephson junctions 1 (IJJs) are of 
great interest nowadays because of its possible appli- 
cations in terahertz physics 2 , quantum electronics, etc. 
Macroscopic quantum tunneling (MQT) in high-T c JJs, 
in particular in IJJ stacks of E^S^CaCuOs+a (BSCCO), 
has been reported in recent years 3-10 , and has attracted 
much attention 11-16 . An unconventional quantum phe- 
nomenon, an enormous enhancement of the MQT rate 5 
due to interaction among stacked IJJs, has also at- 
tracted considerable attention 13-16 . MQT has been stud- 
ied not only in IJJs, but also in various systems, such as 
long asymmetric Josephson junctions and buckling nano- 
bars 17-20 . 

Let T* denote the crossover temperature distinguish- 
ing the MQT and the thermal activation regimes. The 
high T* of IJJs (due to their high critical current density) 
makes them quite attractive as a candidate for a quantum 
bit (qubit) that can be operated at higher temperatures. 
From this point of view, IJJ stacks of La2-^Sr 2 ,Cu04 
(LSCO) compounds could be even more preferable than 
BSCCO stacks, since the Josephson plasma frequency in 
LSCO materials (and, consequently T*) ranges from a 
few hundreds GHz to ~ 1 THz, 21 which is higher than the 
typical value ~ 100 GHz in BSCCO. Moreover, LSCO IJJ 
stacks are also an interesting system to investigate the 
enhancement of the MQT rate 5 , because its inter-layer 
distance (s = 0.7 nm) is shorter than that of BSCCO 
{s = 1.5 nm), so that the interaction among the stacked 
IJJs is expected to be stronger. 

Although the intrinsic Josephson effect in LSCO ma- 
terials has been studied in several papers 22-24 , in most of 
them the sample sizes were not small enough for reliable 



MQT studies. This was due to the technical difficulty to 
fabricate small, submicron, LSCO IJJ stacks carved out 
of a bulk single crystal. Thanks to the recent progress 
achieved in Ref. 25, submicron single-crystal LSCO IJJ 
stacks can now be obtained. 

In this paper, we report precise measurements of 
switching currents in a submicron LSCO IJJ stack in 
which uniform switchings 5 occur. The switching cur- 
rent distributions were measured in a wide temperature 
range, and the crossover temperature T* was determined 
to be T* ps 2 K. In the thermal activation regime, strong 
phase diffusion (re-trapping) has been observed because 
in LSCO the quasi-particle conductivity is high enough 
to suppress the escape rate. We use the theory of MQT in 
IJJs developed in Refs. 14-16, which successfully explains 
the large enhancement of MQT escape rate observed 5 in 
BSCCO IJJ stacks. Due to the interaction among super- 
conducting Cu02 layers in the LSCO IJJ stack, it turns 
out that the phase difference is spatially inhomogeneous, 
both in the inter-plane and the intra-plane directions of 
the IJJ stack. This leads to the enhancement of the es- 
caper rate compared to the standard theory of MQT in 
Josephson junctions. We generalize this approach taking 
into account dissipation and the effect of re-trapping of 
the phase difference. 



II. DEVICE AND SETUP 

A submicron LSCO IJJ stack was fabricated us- 
ing the three dimensional focused-ion-beam etching 
technique 23-25 in a bulk Lai.giSro.ogCu04 single crystal 
grown by the traveling solvent floating zone method 26 . A 
schematic of the LSCO IJJ stack is depicted in Fig. 1(a). 
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The lateral dimensions of the stack are 0.45 x 0.95 /im 2 , 
and the thickness is about 35 ran (with an estimated num- 
ber of layers N « 50). 25 The IJJ stack was placed in a 
copper box inside a 3 He cryostat with a base tempera- 
ture 0.4 K. DC transport measurements were carried out 
with a four-terminal configuration. The low-frequency 
RF noise was filtered by a series of low-pass filters, and 
the high-frequency (~ GHz) thermal noise was attenu- 
ated by a lossy coaxial cable placed between the 1 K pot 
and the sample box. 

Switching currents of the LSCO IJJ stack were mea- 
sured using a high-resolution ramp-timc-based setup, 
which is described in detail elsewhere 7,10 ' 27 . The bias 
current was linearly ramped up at a constant rate 
9.17mA/s, with a repetition rate 21 Hz. In our setup, 
similar to Refs. 10,27, the ramped bias current is turned 
off immediately after detecting a voltage signal above 
20 jitV from the IJJ stack, in order to minimize self- 
heating of the junctions. Switching current distributions 
P(I) have been constructed from 10, 000 switching events 
at each temperature ranging between 0.4 K and 10 K. 

Here we note that the P(I) of a BSCCO IJJ stack 
having a critical current two orders of magnitude smaller 
(~ fxA) than that of the LSCO IJJ stack (~ 100 fj,A) 
was also measured using the same setup, and the results 
(not shown here) were in excellent agreement with the 
thermal activation theory for conventional single Joseph- 
son junctions. In addition, we confirmed that the self 
heating is negligible at 21 Hz, the repetition rate used in 
this measurement, by comparing P(/)'s taken at differ- 
ent repetition rates. These facts imply that a reliable set 
of data can be obtained in our measurement setup. 



III. SWITCHING DYNAMICS OF LSCO IJJ 
STACK 

A current- voltage (/-V) characteristic of the LSCO IJJ 
stack is shown in Fig. 1(b). Unlike typical I-V curves 
of BSCCO, 1 ' 3 ' 5 ' 7 ' 8,10 no clear multi-branches appear, but 
all the junctions simultaneously switch to their voltage 
states. This simultaneous switching event is because all 
the junctions in the LSCO IJJ stack have the same crit- 
ical current, similar to the "uniformly switching" events 
reported by Jin et al. 

The results of the switching measurement are shown 
in Fig. 2(a). The distribution P(I) becomes wider when 
increasing the position of the peak, as the temperature 
decreases from 10 K down to ~ 3.5 K. To see this behav- 
ior of P(I) more clearly, the width (standard deviation) 
(7 of P(I) is plotted as a function of temperature T in 
Fig. 2(b). As shown in the regime (III), <x increases with 
decreasing temperature: from ~ 10 K down to ~ 3.5 K. 
Here we note that the mean switching current < I sw > 
of the LSCO IJJ stack obeys the Ambcgarkar-Baratoff 
formula as shown in Fig. 2(c), indicating that our LSCO 
IJJ stack consists of an array of tunnel junctions. 

From 3.5 K down to 2K, the width a of P(I) shrinks 
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FIG. 1: (Color online) (a) Schematic of the geometry of the 
sample (left) and an array of LSCO IJJs (right), (b) A 
current-voltage (I-V) characteristic of the IJJ stack measured 
at 4.2 K. All junctions in the stack have the same critical cur- 
rent, and therefore simultaneously switch to the voltage state 
(arrows) without showing multi-branches. 



very slightly, as seen in the temperature regime (II) in 
Fig. 2(b). Finally P(I) does not change below ~ 2 K, and 
therefore a saturates, as seen in (I) in Fig. 2(b). The sat- 
uration of a in the lower-temperature regime (I) qualita- 
tively suggests that MQT occurs in these temperatures. 
The temperature dependence of a in the regime (III), 
however, does not correspond to the power law a oc T 2 / 3 , 
which is expected in the thermal activation theory of con- 
ventional single junctions 27 . Thus the switching dynam- 
ics of the LSCO IJJ stack cannot be understood only 
by the conventional MQT and thermal activation the- 
ories for single junctions. In the following sections we 
quantitatively discuss those switchings' dynamics, which 
are in good agreement with the MQT theory for stacked 
IJJs 14-16 and thermal activation theory in the presence 
of moderate damping 28,34 . 



IV. THEORETICAL MODEL AND ANALYSIS 

Since the position of the peaks in the current distri- 
bution P(I) does not change at low temperatures [Fig. 
2(a)], we assume that the switching events are dominated 
by MQT. The crossover temperature T* between MQT 
and thermal-activated regimes is approximately 2 K. For 
temperatures larger than approximately 3.5 K, a starts 
to decrease when the temperature increases. Let us first 
focus in the low-temperature regime T < T* . Afterwards 
we will discuss the thermal-activated regime, especially 
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FIG. 2: (Color online) (a) Switching current distributions 
P(I) at different temperatures, (b) The width (standard de- 
viation a) of P(I) as a function of temperature T. The ver- 
tical dashed lines qualitatively separate the different temper- 
ature regimes: (I) MQT, (II) cross-over, and (III) thermal 
activation with phase re-trapping (phase diffusion), (c) Mean 
switching current < 7 SW > versus T. The dots are exper- 
imental data, and the dashed curve is a calculation of the 
Ambegaokar-Baratoff theory 31 with assuming that the super- 
conducting gap is BCS type. 



the unusual behavior of a with temperature, by taking 
into account the phase re-trapping (phase diffusion) ef- 
fect due to dissipation. 



A. Regime of macroscopic quantum tunneling 

In the MQT regime, it is more convenient to consider 
the escape rate instead of the switching probability 
P(I). The escape rate T(I) can be calculated from P(I) 



T(I) 



(dl/dt) P(I) 

l-fdPP(F] 
o 



(1) 



where dl/dt (= 9.17 mA/s) is the rate of current ramp 
in the switching measurements. The data for T(I) cal- 
culated from P(I), measured at T = 0.4 K, is shown in 
Fig. 4. The standard theory of MQT in Josephson junc- 
tions gives the following formula for the escape rate 29 : 



Tmqt = ^aVUOirB [1 - (I/I c ) 2 ] 1/4 exp(-B), (2) 

where I is a bias current, I c is the critical current of the 
junction, uj p is the plasma frequency at zero bias cur- 
rent, and the prefactor a is of the order of unity (for 
more details about a, see, e.g., the review 30). Also 
B = S/h, where S is the action of the system calcu- 
lated for an imaginary-time trajectory of the phase dif- 
ference <p(r) (bounce solution). For a single Josephson 
junction in the short junction regime, i.e., Xj > D, where 
Aj the Josephson penetration depth and D the junction 
length, the phase difference <p(t) can be considered as a 
coordinate of a particle, moving in the inverse potential 
— U (p) oc I c cos <p + Ip. For a current / close to a critical 
current I c , and in the absence of dissipation, the result 
for the tunneling exponent is 29 
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12J C [1 - (I/I c ) 2 ] 
5eu p (I/I c ) 2 
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(3) 



For the LSCO IJJ stack under study, I c = 344 fj,A and 
f p = ujp/2-K = 440 GHz (sec Appendix). Substituting 
these values into Eqs. (2) and (3), we obtain the result for 
T(I) which is many orders of magnitude smaller than the 
experimental value for all currents measured. In other 
words, the formula (3) strongly overestimates the real 
value of the tunneling exponent B. Thus, we conclude, 
that the sample under study cannot be considered as a 
single Josephson junction. 

To explain the unusual large escape rate T both in the 
MQT and the thermal-activated regimes, we will exploit 
the theory of MQT developed in Rcfs. 14-16. Namely, we 
model the sample as a stack of intrinsic Josephson junc- 
tions formed by Cu02 superconducting layers, and tp n be- 
ing the phase difference between nth and (n + l)th layers 
in the stack. The number of stacked IJJs was estimated 
to be N « 50 by measuring the thickness of the IJJ stack 
and the voltage gap in an I—V curve. The schematic of 
the sample is shown in Fig. 3. A DC-bias current flows 
along the y-direction (perpendicular to Cu02 planes) in- 
side the IJJ stack, whereas it flows along the x-direction 
in all other places (see arrows in Fig. 3). The j/-axis is 
perpendicular to the area of the IJJ stack. The phase dif- 
ferences (p n are assumed to be spatially inhomogeneous 
in the x-direction. 

We represent <p n in the form ip n = po + ip n (x,t), where 
Po = arcsin(/// c ) corresponds to the equilibrium value 



4 



o 

e 



i/j (x,t) 



Xq 



535 



yo 



D 



X 



FIG. 3: (Color online) Schematic of the sample considered 
in the theoretical model. The phase difference ip(x, y, r) is 
assumed to be ^-independent in some region in the middle of 
the stack, and exponentially decreases outside it. The phase 
difference ipo(x, r) is inhomogeneous in the s-direction, giving 
the characteristic size xo of the tunneling fluxon. 



of the phase difference. The tfj n satisfy the equation of 
motion (for a review, see, e.g., Ref. 2) 
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dt 
dx 2 



1) 



0. 



(4) 



where j = I/I c is a normalized current, s is the inter- 
layer distance, A a h is the in-plane penetration depth, 
a = 4:Tra±/(eujp) is the dissipation parameter, where 
e is the dielectric constant of the insulating layer, and 
ct_l is the quasiparticle conductivity across the layers. 
In Eq. (4), the time t is normalized by l/(J p , the in- 
coordinate is normalized by the out-of-plane penetra- 
tion depth A c , and the action of the d 2 operator is 
d 2 f n = fn+i — 2/„ + fn-i- We will use below the contin- 
uum limit, replacing ?p n (x,t) by the continuum function 
ij)(x,y,t), where the y-coordinate (perpendicular to the 
layers) is normalized by X a {,. In this case, we can replace 
iKJs 2 )d 2 n d 2 /dy 2 . 

To estimate the dissipation parameter a = 
4:TT(j±/(eujp), we use the Ambegaokar-Baratoff for- 
mula for the critical current 31 (which works well for 
LSCO 32 ), given by I c s/(cr ± A) = 7rA /(2e), where A 
is the junction's area and A is the superconducting 
gap at zero temperature. Using the BCS relation 



Ao = 1.76fcsT c , where ks is the Boltzmann constant 
and T c is the superconducting transition temperature, 
and taking the value of the dielectric constant 21 e = 25, 
we obtain the estimation a ~ 0.2 corresponding to 
moderate damping (cf. a -C 1 for BSCCO). 

To find the escape rate T one should find a periodic 
solution to the equation of motion for ij)(x, y, r) in imag- 
inary time t = vt, with period K/ksT, and calculate 
the action S of the system corresponding to this equa- 
tion (which is actually non-local in both the y-coordinatc 
and the imaginary time r). This is a hard numerical 
task. Here, we use an approximate variational approach 
similar to that developed in Ref. 16. There, it is as- 
sumed that the tunneling occurs mainly in one junction 
in the stack, say at point yo, creating exponentially de- 
caying tails of the phase difference outside it, that is 
ip(x,y,r) oc exp(— n\y — y \) (see Fig. 3). The analysis of 
the linearized Eq. (4) shows that n ss n\ c /[D(l — j 2 ) 1 / 4 }, 
where D is the junction width in the x-dircction. For 
the BSCCO samples considered in Ref. 16, Ks/A a f, > 1. 
This means that the phase difference indeed decreases 
fast with increasing y, and one can consider the tunneling 
through one junction. By contrast, we have Ks/X a b 1 
for the LSCO IJJ stack under study. It actually reflects 
the fact that in LSCO s is smaller than that of BSCCO, 
so that the interaction among Cu02 is stronger. 

Here, we assume the following profile for the phase dif- 
ference ip(x, y, t) in the {/-direction: %[){x, y, r) = ipo(x, r) 
inside the region \y — y a \ < n = n (l — j 2 ) 1 / 4 , and 
ip(x,y,r) exponentially decays outside this region (see 
Fig. 3). The parameter no ~ D\ a b/ (2-k\ c s) can be con- 
sidered as the effective number of junctions taking part 
in the tunneling process. Taking the following estimates: 
the anisotropy parameter 7 = \ C /X a b = 30, s = 0.7 nm, 
and D = 0.95 /im, we obtain no pa 5. Note also the ad- 
ditional factor (1 — j 2 ) 1 / 4 , which increases the exponent 
k and makes the number of junctions taking part in the 
tunneling smaller for / closer to I c . 

Now we obtain the effective action for the phase dif- 
ference ipo( x i T ) m a way similar to that used in Ref. 16. 
We solve the linearized Eq. (4) at |y — j/o| > ^o, with the 
boundary condition Tp(yo,x,r) = -0o(x,t). Using then 
Eq. (4) and Maxwell equations, we obtain the relation 
for i/jq(x,t) and the action corresponding to it (for de- 
tails, see Ref. 16). As a result, we find: 

Seff = ^ Jdrjdxl^ (^-j +Mo(l - cosi/>o) 

^ 

d 

-j(^o - sinV>o) + ^7^T~ir(dx'P{x,x')^^- 



4ttD dx 



a_ r , [j} {x,T) -i/j (x,t')Y 
4tt J T (t-t>) 2 

— OO 



(5) 



where (Xq = y/l — j 2 , d = D/X c , tq = TUOp/ksT, A 
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hI c no/(2eojp), and 



P(x,x') = In 



sin[7r(x + x')/d] 



sin[7r(a; — x')/d] 



(6) 



The action (5) is similar to the action of a single (long) 
Josephson junction 15,33 , but it contains a term non-local 
in the x-coordinate, instead of a local term proportional 
to (dip/dx) 2 . This non-locality, originating from the in- 
teraction among many junctions in the stack, strongly 
reduces the (positive) contribution to the action from 
the ^-derivatives of the phase difference, compared to 
the single- Josephson junction. This favors the phase dif- 
ference ipo to be spatial inhomogeneous for stacks of the 
width D exceeding several 7s (instead of Aj for a single 
JJ, Ref. 15). 

In further approximation, we represent tPq(x,t) in the 
form ipo(x,r) = F(x)q(r), where F(x) is a trial func- 
tion, describing the shape of the tunneling fluxon in the 
x-direction, and q(r) plays a role of the collective coor- 
dinate. The form of the function F(x) is chosen from 
the physical reason that the fluxon nucleates at the junc- 
tion's edge 14-16 . We considered different trial functions, 
and obtained similar results. For all results in this paper, 
the function F(x) used is 



F(x) 



C(x Q ) 



1 + (x/x y 



(7) 



where Xq is the characteristic size of the tunneling fluxon 
(see Fig. 3), and the normalization constant C(xq) is cho- 
sen such that Jq dxF 2 (x)/d = 1 . The analysis shows that 
xq is about several 7s (for more details, see Ref. 16). 
Since for the LSCO IJJ stack under study xq ~ 7s sa 
20 nm -C D = 0.95 /Ltm, the inhomogeneity of the phase 
difference in the x-dircction is crucial to make a large 
difference from the conventional MQT theory. 

Substituting Tpo(x,r) in the form tpo(x,T) = F(x)q(r) 
into Eq. (5) and expanding the expression inside the in- 
tegrals in a Taylor series, we obtain a particle-like action 
in the form: 



10 9 



No dissipation 




• 0.4 K 



305 



310 315 
Current / (uA) 



320 



FIG. 4: (Color online) The escape rate V versus current I. 
Dots correspond to the experimental data at T = 0.4 K. The 
solid curve is the theoretical prediction. The parameters of 
the model are: 7 = 30, no = 5, s = 0.7 nm, D = 0.95 /im, 
J c = 344 uA, f p = 440 GHz, T = 0, and a = 0.3. The 
dashed curve is calculated for zero dissipation, a = 0. Other 
parameters are the same as for the solid curve. 



and 



dxF n {x). 



(11) 



The analysis shows that for the parameters under study, 
we can neglect in Eq. (9) the terms with q e and higher. 
The potential U (q) implicitly depends on the fitting pa- 
rameter xq. The collective coordinate q(r) in Eq. (8) 
satisfies the imaginary-time equation of motion 



dU 
dq 



dT 



, g(r)-g(r') 
(r - r') 2 



(12) 



T OO 



where the potential U(q) is 



U(q) 



Mo 



E 

71=2 



(2n)! q 



(-l) n l/ 2n +l 2n+l 

{ (2,1 + 1)! q 



In this equation, 
jj, n s 



d d 



/! = /J + 



47rA„, ) tf 2 



dx 



(8) 



(9) 



with periodic boundary conditions q(r + tq) = q(r). We 
solve this equation numerically for given xo, and calcu- 
late the tunneling exponent B = S e g/H. We repeat this 
procedure until finding the xq corresponding to the min- 
imum of B. The MQT escape rate is calculated then by 
Eq. (2). 

At temperatures much smaller than the crossover tem- 
perature T* , the escape rate T is independent of T. Fig- 
ure 4 shows the result of the calculations of T(I) at 
T <C T* with the prcfactor a = 1. Here the dissipa- 
tion parameter a = 0.3, which is in reasonable agree- 
ment with the estimated value 0.2, is found to be the 
best fit. Note the strong dependence of the tunneling 
exponent B, and consequently of the MQT rate T on the 
dissipation parameter a. Even a not too high dissipation 
constant a = 0.3 strongly reduces the escape rate T (see 
the dashed and the solid curve in Fig. 4). 
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FIG. 5: (Color online) Switching current distributions P(I) at 
three different temperatures above T*. Dots are experimen- 
tal data. Dashed curves are calculated according to Eq. (16). 
The parameters of the model are the same as for Fig. 4. 
Solid curves are calculated taking into account the phase re- 
trapping effect. 



B. Thermal-activated regime 

Let us now focus on the thermal-activated regime. At 
temperatures above the crossover temperature T*, the 
solution of Eq. (12) with periodic boundary conditions 
q(r + tq) = q(r) becomes trivial: q(r) = qo, where qo 
corresponds to the maximum of the potential U(q). The 
crossover temperature T* is denned as 



hujp 

2nki 



TT „, a a 



(13) 



where U is the second derivative of U (q) at q — qo . The 
theoretical prediction of T* turns out to be ~ 2K for 
the bias currents / used in the measurement with the 
same parameters as in Fig. 4, even though T* depends 
on the bias current /. This prediction of T* is in good 
agreement with the experiments [see Fig. 2(b)]. 
For T > T*, the tunneling exponent B reduces to 

AujpUo 



B = 



(14) 



where Uq = £7(<7o) is the maximum of the potential U(q). 
To calculate the thermally activated escape rate, the 
value of Uq has to be minimized with respect to the char- 
acteristic size of the fluxon xq. The tunneling exponent 
B in the thermal-activated regime does not depend on 
the dissipation a, and the escape rate (2) only slightly 
depends on a due to the prefactor a. For T > T* , the 
latter one can be written as 29,30 , 



JL 

Wo 



n 



au) r . 



/' 



auj r , 



(15) 



where ui n — 2im/TQ. 

Now it is more convenient to consider the switch- 
ing current distribution PtaCO in the thermal-activated 
regime instead of Tta(-0- Using Eq. (1), we obtain 



Pta(I) = 



Tta(/) 
(dl/dt) 



exp 



(dl/dt) 



Jdl'TTA(l') 



(16) 

The results of our calculations of Pta(E) at three different 
temperatures T > T* are shown by the dashed lines in 
Fig. 5. For all curves, the parameters are the same as for 
Fig. 4. We see good agreement with experimental data 
at T — 4K. For larger temperatures, however, the peaks 
of the theoretical curves shift to lower currents, in com- 
parison to the experimental curves. In addition to that, 
the width of the theoretical Pta(-0 increases as the tem- 
perature increases, while the experimental _Pta(-0 shows 
the opposite behavior. We attribute these small discrep- 
ancies to the phase re-trapping (phase diffusion) process 
due to dissipation. In the next section, we will address 
and solve this issue. 



C. Thermal activation with phase re-trapping 

Thermal fluctuations stimulate the switching of the 
phase ip of the junction to the running (resistive) state. 
However, similar fluctuations can also help re-trap ip back 
to the metastablc (superconducting) state. For more in- 
formation about re-trapping in Josephson junctions, see, 
e.g., Ref. 28. 

In contrast to the normal switching probability with- 
out dissipation in the thermal-activated regime T > T*, 
the probability of re-trapping strongly depends on the 
dissipation parameter a. We consider the re-trapping ef- 
fect following the theory developed in Ref. 34. There, a 
dissipative particle, moving in a periodic washboard po- 
tential is studied in the presence of thermal noise. The 
state of the particle with a constant average velocity cor- 
responds to the resistive state of the Josephson junction. 
Thermal fluctuations can trap the particle in one of the 
local minima of the potential (corresponding to the su- 
perconducting state). The authors of Ref. 34 derived an 
analytical formula for the re-trapping rate r r . Applying 
this to the Josephson junction, T r is then given by 



' Ej(i-i r y 

27rr 2 fc B T 



■ exp 



Ej(I - Ir? 

2I?a 2 k B T 



(17) 



where I r = Aal c /ir is the re-trapping current, and Ej = 
hl c /(2e). This formula is valid for a single short Joseph- 
son junction. In our case the situation is much more com- 
plicated because the phase difference is inhomogeneous in 
the a;-direction. We assume that the phase difference can 
be "partially re-trapped" . In this process, thermal fluc- 
tuations lead to the appearance of a "rc-trapped fluxon" 
of width x r in the cc-direction. After this, the junction 
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dynamically relaxes to the mctastablc (superconducting) 
state. We modify Eq. (17) by multiplying the Josephson 
energy Ej by the factor x r /D, where we take x r equal 

to Xq. 

Using Eq. (17) for the re-trapping rate (with an addi- 
tional factor xq/D in Ej), we calculate the probability 
Fur{I) for the fluxon not being re-trapped. The rela- 
tionship between T r (I) and F n n(I) becomes 28 : 

Fn R (I)=e^^-^jdI'T r (I^ . (18) 

Thus, the switching distribution, which is actually mea- 
sured in the experiment, is equal to Pta(-0-Fti.r(-0, where 
-Pta(-0 is given by Eq. (16). The results of such calcula- 
tions at three different temperatures, from 4 to 6 K, are 
shown in Fig. 5 by the solid curves. For all curves, the 
model parameters are the same as for Fig. 4. Now we see 
that the theoretical curves show better agreement with 
the experimental data: the positions of P{I) do not shift 
towards smaller currents, and the widths of the P(I) de- 
crease with increasing temperature. For T = 4K, the 
solid and dashed (no re-trapping) curves practically co- 
incide. This means that the re-trapping effect turns out 
to be not significant for T < 4 K. This correlates well 
with the temperature dependence of the standard devia- 
tion a [width of P(I)] in Fig. 2(b). Namely, due to the 
existence of phase diffusion at T > 4 K, <r(T) shows un- 
conventional temperature dependence, different from the 
T 2 / 3 law. This phenomenon [violation of the T 2 / 3 law of 
ct(T)], unique for the LSCO IJJ stack [cf. Refs. 3,5,7,10 
(BSCCO) and Ref. 9 (HBCCO)], can be attributed to the 
relatively small superconducting gap of LSCO; namely, 
above 4 K there are still non-negligible quasi-particlcs 
which contribute dissipation to the IJJ stack, resulting 
in re-trapping. 



V. CONCLUSIONS 



Acknowledgements 

The authors thank S. Kawabata, A. Tanaka, 
T. Koyama, M. Machida, I. Kakeya, N. Kobayashi, and 
Y. Ootuka for discussions. This work was partially sup- 
ported by MEXT under Grant No. 1905014 and the 
Russian Foundation for Basic Research (projects Nos. 
09-02-92114 and 10-02-92600). Y.K. was financially sup- 
ported by NIMS. A.O.S. acknowledges support from the 
Dynasty Foundation. F.N. was partially supported by 
LPS, NSA, ARO, NSF grant No. 0726909, Grant-in- Aid 
for Scientific Research (S), MEXT Kakenhi on Quantum 
Cybernetics, and the JSPS via its FIRST program. 




300 310 320 330 
Current / (^A) 



FIG. 6: (Color online) Determination of I c from the experi- 
ments. Dots are experimental data converted using Eq. (A2), 
and dashed lines are linear fits. The mean value of I at the 
intercepts of the extrapolated lines gives I c — 344 fiA. 



We have measured the switching current distributions 
in a submicron LSCO IJJ stack for a wide temperature 
range, from 0.4 K to 10 K. The switching probability 
does not depend on temperature for T < 2K, implying 
that the LSCO IJJ stack is in the MQT regime. We 
have applied and extended the theory for switching dy- 
namics in IJJ stacks taking into account the effect of 
dissipation both in the MQT and the thermal activated 
regimes. The theory and the experiment are in good 
agreement; therefore the enhancement of the MQT rate 
due to interactions among IJJs is confirmed to be a uni- 
versal phenomenon in IJJ stacks. We have also shown 
that dissipation plays an important role for all tempera- 
tures. The phase re-trapping takes place in the thermal- 
activated regime, where T > 4K. 



Appendix A: Determination of the fluctuation-free 
critical current 

Here we describe how the fluctuation-free critical cur- 
rent I c was determined from the experiments. Here we 
used the escape rate in the thermal-activated regime in 
the form 

r Ti = g (w¥ exp(-||), (A!, 

where the prefactor a is assumed to be unity. Then 
Eq. (Al) is converted 2 to a linear function with respect 
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to the bias current /, 



-In 



2tiT 



TA 



wp(i-r»)V* 




(A2) 

where the barrier height AU was approximated by a cu- 
bic function 35 : AU « ^£j(l - I/I c ) 3/2 - The normal- 
ized r after this conversion is shown in Fig. 6. Here we 
used the following values as the initial guesses for the 
fit: I c = 340 fiA, obtained by a rough fit with Eq. (Al) 
in which I c and T arc defined as free parameters, and 
C = 1 35 fF determined from the geometry of the junc- 
tion, i.e., C = EqsA/s. 

From the mean value of the zero-crossing of the extrap- 



olated lines in Fig. 6, I c was determined to be 344 fxA. 
Since u) p is inside the logarithmic term on the left hand 
side of Eq. (A2), the variation of the initial I c gives a 
very small difference to the result of the conversion 27 in 
Eq. (A2). Therefore it is sufficient to use the roughly 
estimated I c as an initial value, and consequently I c can 
be determined by repeating the fitting procedure with 
high accuracy. The zero-bias plasma frequency was also 
obtained: u p = v / 2eI c /hC = 2tt x 440 GHz. In Fig. 6, 
the data obtained at temperatures T < 3.5 K have been 
chosen for the fit, because the re-trapping effect cannot 
be negligible at temperatures T > 4K, as discussed in 
Sect. IV C and shown in Fig. 5. 
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